Vortex configurations and metastability in mesoscopic superconductors 
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The vortex dynamics in mesoscopic superconducting cylinders with rectangular cross section under an axially 
applied magnetic field is investigated in the multivortex London regime. The rectangles considered range from a 
square up to an infinite slab. The flux distribution and total flux carried by a vortex placed in an arbitrary position 
of the sample is calculated analytically by assuming Clem's solution for the vortex core. The Bean-Livingston 
energy barrier is also analytically calculated in this framework. A Langevin algorithm simulates the flux penetra- 
tion and dynamical evolution of the vortices as the external field is slowly cycled. The simulated magnetization 
process is governed by metastable states. The magnetization curves are hysteretic, with paramagnetic response 
in part of the downward branch, and present a series of peaks corresponding to the entry or expulsion of a single 
vortex. For elongated rectangles, the vortices arrange themselves into parallel vortex chains and an additional 
modulation of the magnetization, corresponding to creation or destruction of a vortex chain, comes out. 
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1. Introduction 

Since the pioneering work of Abrikosov pp, it 
is well known that in a macroscopic type-II su- 
perconductor the magnetic field, in the so-called 
mixed state, penetrates in the form of singly 
quantized flux lines (or vortices), which form a 
triangular lattice. For finite samples the vortices 
are formed at the surface and then pulled to the 
sample center by the shielding Meissner currents. 
Notwithstanding, close to the sample edge, a vor- 
tex is strongly attracted by the superconductor- 
vacuum interface. These competing interactions 
give rise to a energy barrier, the so-called Bean- 
Livingston surface barrier, which retards the 
movement of vortices towards the sample cen- 
ter |2I3| . In samples with very smooth interfaces, 
this barrier leads to irreversible field-dependent 
magnetization loops M[H) and finite critical cur- 
rents even in the absence of inhomogeneities |4I5| . 
Nevertheless, defects in the interface may cause 
the local destruction of the surface barrier, open- 
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ing leaks for vortex entry. 

Recent progress in nanostructuring of super- 
conducting materials |fil7| has provided the op- 
portunity of studying the vortex state in meso- 
scopic samples, i.e. whose dimensions are of the 
order of the penetration depth A and/or the co- 
herence length £, with sharp interfaces and no no- 
ticeable pinning. At these length scales the vor- 
tex lattice and the vortex itself may present new 
and very interesting properties, as for example, 
the formation of multi-quanta giant vortices, vor- 
tex molecules and chain-like vortex arrangements. 
These vortex structures are strongly dependent 
on the sample geometry and size, as has been 
shown by numerical and analytical solutions of 
the Ginzburg-Landau (GL) equations and 
experiments |6I10) . and may appear even in sam- 
ples made of type-I materials (with k — A/£ < 
1 / y2) , such as aluminium. The reason is that in 
thin films the parameter governing the appear- 
ance of vortices is actually the effective GL pa- 
rameter, k — A/£, where A = X 2 /d may be much 
greater than A for small sample thicknesses d. 
The choice between the multivortex and the gi- 
ant vortex states depends on the value of k |Bj and 
on the system size as well. Here we shall consider 
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only systems with high k and sizes comparable to 
A, though larger than £. For these systems the 
multivortex state persists over a wide area of the 
phase diagram. 

In the course of our recent research with meso- 
scopic superconductors |14I15I16| we have been 
studying the vortex dynamics, within the London 
approximation, in films, multilayers, and strips, 
submitted to an external magnetic field and dc 
currents. In this article we present new results 
for mesoscopic rectangles in the presence of an 
applied magnetic field. The vortex structure and 
the energy surface barrier present in these sys- 
tems are described within the London limit of the 
Ginzburg-Landau equations. The magnetization 
process is simulated by a fast algorithm based on 
Langevin dynamics. We observe, for elongated 
rectangles, the appearance of vortex chain states 
and transitions between these states resulting in a 
modulation of the magnetization loop. Particular 
attention is given to the formation of metastable 
vortex configurations in these systems during the 
magnetization process, which leads to hysteretic 
magnetization curves with paramagnetic response 
in part of the downward branch. 

2. Vortex structure and the Bean- 
Livingston barrier 

In the high-K limit, the superconducting order 
parameter is essentially homogeneous except near 
a vortex core, where its spatial distribution may 
be given approximately by |1 lj 
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where p is the distance to the vortex center. For 
points r = (x, y) far away from the vortex core, 
located at r' = (x',y'), \ip(r)\ is uniform and the 
second Ginzburg-Landau (GL) equation reduces 
to the London equation , 
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The solution near the vortex core may be accom- 
plished very precisely by inserting the variational 
trial function into the GL free energy func- 
tional. The result is equivalent to making the 



cutoff b(r) — > b(y / r 2 + 2£ 2 ) in the London solu- 
tion for the vortex local induction. 

In this spirit, one may compute the flux dis- 
tribution of a vortex confined in a cylinder if the 
appropriate boundary condition is used. Here we 
consider long cylinders with a rectangular cross 
section of width W and length L. An exter- 
nal magnetic field H is applied axially and the 
vortices are assumed to be perfectly aligned with 
H. This problem has been considered previously 
within the London theory by Sardella, Doria and 
Netto 521- Here we shall use Clem's variational 
solution for the vortex core to compute the local 
field distribution and the position-dependent, ef- 
fective flux carried by a vortex. The result for the 
local induction generated by a vortex is 
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where k m = rmr/W, a + = a + a' and a_ = 
| a — a'\ + £ (a — x,y). With this solu- 
tion, one can compute the interaction energy be- 
tween vortex i and j in the usual way, that is, 
E vv (xi ,yi',Xj,yj) — fob^ {x t ,y t ). Integrating the 
local induction b v {x,y) over the entire sample 
area, one finds the effective magnetic flux carried 
by a vortex, 

cosh(:r/A - W/2X) 4 



(x,y) = 0o [l- 



sin k. 



coshW72A 
% x cosh k m (y 



X 2 W 



m—l 



" m 



cosh k m L/2 



(4) 



The prime in the sum operator indicates that 
only the terms where m is odd are taken into 
account. Note that in mesoscopic samples the 
effective flux may be quite smaller than the flux 
quantum, even for a vortex positioned at the sam- 
ple center, where 4>(x, y) is maximum. 

The homogeneous solution of Eq. [^corresponds 
to the local Meissner screening flux distribution, 
which is given by 
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This field distribution is minimum at the sam- 
ples center. Thus, the Meissner screening current 
density, jm = /i 'Vx 1>m, repels the vortex away 
from the surfaces, towards the sample center. The 
corresponding potential energy felt by the vortex 
is then given by E M = (f> b M (x', y'). 

On the other hand, there is an energy cost to 
put a vortex inside the superconducting speci- 
men. This is given by the vortex self-energy 
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This energy depends on the vortex position and is 
maximum at the sample center. That is, the vor- 
tex is attracted by the surfaces or, in other words, 
by its images, which are necessary to satisfy the 
boundary conditions. The attractive self-energy 
and the repulsive Meissner energy form the well 
known surface (Bean-Livingston) barrier, which 
delays vortex entrance and exit in the supercon- 
ducting sample. 

In Fig. 2] we make a surface plot of the en- 
ergy distribution in a superconducting rectangle 
with W = 0.2A = 20£ and L = 2W for three 
different situations. In Fig. ^a), H equals the 
lower critical field, H = H c i, for which there is 
a global minimum at the sample center and a 
strong energy barrier near the sample surfaces. 
In mesoscopic samples the activation energy to 
overcome this barrier is usually much higher than 
the thermal energy ^T. In Fig.^b), H = H en , 
where the surface barrier for the first vortex en- 
try disappears. The entrance field may be de- 
fined as dE se if/dx = —dEM(H en )/dx, with the 
derivatives being evaluated at one of the points 
(x, y) = (£, L/2) or (D - £, L/2), as suggested by 
Fig. ^b) . Note that in the smaller edges of the 
sample, the barrier is still present, which means 
that in a rectangle the first vortex tends to enter 
the sample through one of the larger edges. The 
total energy when a vortex is placed at the sam- 
ple center at H = H en is depicted in Fig. ^c) . In 
this case, a new energy barrier for the entrance 



of a second vortex is developed. This means 
that, at zero temperature, the second vortex will 
be allowed to enter the sample only at a field 
H2 > H en . In the field region H en < H < H2, 
the sample behaves as in the Meissner state. 
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Figure 1. (a) and (b): Energy distribution in a 
superconducting rectangle with W = 0.2A = 20£ 
and L = 2W in the absence of vortices at the first 
critical field, H = H c \ (a), and at the entrance 
field, H = H en (b). (c) The total energy when 
a vortex is located at the sample center at H = 

H en . 
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3. Langevin dynamics and simulation 
scheme 

Next, we review our model to study vortex dy- 
namics in confined geometries. The time evo- 
lution of a vortex is described by overdamped 
Langevin equations of motion, 

V v i = -V i E + Y i {t), (7) 

where 

Y,[Eself{ri)+E M {ri)] (8) 

i 

is the total energy of the vortex distribution, r\ is 
the Bardeen-Stephen friction coefficient, Vj is the 
vortex i velocity. is a Gaussian stochastic noise 
related with a small temperature T and the fric- 
tionTyby (T ati (t)Tp d (t')) = 2r)k B T5 a f } S ij 5(t-t'), 
where the Greek indices stand for the directions 
x or y. This temperature plays the role of ac- 
celerating the convergence towards a stationary 
state. 

The simulation consists in numerically inte- 
grating the coupled Langevin equations for all 
vortices inside the sample. The integration is 
based on a finite difference algorithm. A vor- 
tex is allowed to enter the sample and partici- 
pate on the dynamics if it satisfies a force balance 
condition near one of the sample surfaces. The 
procedure is the same adopted in earlier publica- 
tions [14115116)] : at each time step a test vortex is 
placed in a random position a distance £ from one 
of the sample interfaces and the total force act- 
ing on it is calculated. If this force points to the 
sample interior, the vortex is accepted, otherwise 
it is rejected. 

The magnetization loops are calculated as an 
external field is varied in small steps during up 
to 5 x 10 5 time steps. We verified that the time 
between successive field changes is large enough 
for the system to relax to the desired station- 
ary states. The magnetization M = ^q 1 B — H 
is calculated by integrating the total amount of 
magnetic induction B inside the sample, B = 
^[Jdxdyb M (x,y) + E,: <X^ 2/0] • 



4. Vortex chain states in mesoscopic films 

Here, we consider the case where L — ► oo, 
which corresponds to a film infinitely long in the 
y and z directions. To perform the simulation, we 
take a slice L y in the y direction and assume peri- 
odic boundary conditions at y = 0, L y , in such a 
way that the simulation is restricted to the region 
Q < x < W, < y < L y . 

In Fig. Ufa) we show magnetization loops for 
several film thicknesses, W/X = 1, 2 and 4. We 
choose k — 20 and L y — 40A. All the loops are 
characterized by hysteretic behavior and a distri- 
bution of peaks at specific field values. Snapshots 
of the vortex configurations [some of which are de- 
picted in Fig. Elb)] show that the VL is composed 
of n linear chains of vortices parallel to the film 
surfaces and the extra magnetization peaks are 
associated with sudden rearrangements of the VL 
to accommodate a new vortex chain (n — > n + 1 
transitions), in the case of increasing magnetic 
field, or destroy an existent chain (n — > n — 1 
transitions), for decreasing field. Note that as 
the film thickness increases the matching peaks 
become less evident in the magnetization loops, 
that is, the film gradually crosses over from meso- 
scopic to macroscopic behavior. The matching 
effect of vortex chains in thin slabs (W < A) has 
been studied in the framework of thermodynamic 
equilibrium^]. Our approach assumes that a 
vortex nucleates at the surfaces, i.e., for fields just 
above H c \ it has to overcome a strong energy bar- 
rier in order to find the global minimum at the 
equatorial plane of the film. Hence we are deal- 
ing with steady metastable states. This situation 
is closer to real experimental conditions, where 
the time necessary for the system to relax to the 
thermodynamic equilibrium is not accessible. In 
this case the critical fields where the transitions 
n — > n ± 1 take place are history dependent. In 
Fig. Etb) we show snapshots of the vortex config- 
uration in a slice of the W = A film for different 
points of the magnetization cycle indicated in the 
upper panel of Fig. Ufa). 
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Figure 2. (a) Field dependent magnetization 
loops for infinite films of width W = A, 2A and 
4A. The dashed lines are just guides for a bet- 
ter identification of the transitions between suc- 
cessive vortex chain states, (b) Snapshots of the 
vortex lattice in the points of the magnetization 
loop indicated at the W = A film magnetization 
curve. 
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Figure 3. (a) Magnetization loops of rectangles 
of different aspect ratio, W/L = 1/1, 1/2, 1/4, 
and 1/8, with W = 0.2A = 20£. 



5. Vortex states in mesoscopic rectangles 

Now we consider the more general case of long 
cylinders with rectangular cross section. The cal- 
culations were performed for different aspect ra- 
tios, W/L = 1/1, 1/2, 1/4, and 1/8, with W fixed 
at W — 0.2A and k = 100. The magnetization 
loops for these samples are depicted in Fig. [3] and 
in Fig. 0] we show some vortex configurations for 
the L = W and L — AW samples. The magneti- 
zation peaks correspond to entry (in the upward 
magnetization branch) or expulsion (in the down- 
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ward branch) of an individual vortex. In some 
points of the magnetization curve of the longer 
samples (L = AW and L = 8W), jumps of two 
or more vortices are also observed, specially at 
higher fields. 

States with same vorticity may have two dis- 
tinct metastable solutions, depending on whether 
the field is going up or down. This is illustrated 
in Figs. 0Jc) and (d), where in the state with 
vorticity five in a square, vortices arrange in a 
pentagon, in upward field, or in a face centered 
square, in downward field. 

As the rectangle is elongated the vortices or- 
ganize into parallel chains of vortices resembling 
the chain states in infinite films. The effect of 
these arrangements are seen as a modulation of 
the magnetization loop. Note that the magnetiza- 
tion curve of the L = 8W rectangle is particularly 
similar to that of the magnetization of the W = A 
infinite film. In both systems the ratio W/^ = 20. 

As illustrated in Figs. and hysteresis in 
the magnetization curve is a property shared by 
all the samples studied here and is a result of 
the different circumstances in which vortices en- 
ter or are expelled from the sample through the 
surface barrier. The surface barrier delays the 
incursion of vortices towards the sample interior. 
As a result, in increasing magnetic field, the vor- 
tex lattice is driven into successive superheated 
metastable states and the diamagnetic response 
is stronger than the expected for the thermody- 
namic equilibrium. In decreasing field, a surface 
barrier is also present, now against vortex ex- 
pulsion. In this case the system is driven into 
supercooled metastable states and the response 
is paramagnetic in most part of the downward 
branch of the magnetization loop. It is interest- 
ing to note that this behavior of the magnetic 
response, diamagnetic in upward field and para- 
magnetic in part of the downward branch, is quite 
frequent in superconducting mesoscopic samples 
of different geometries and in a broad range of k 
values |17I8I18| . It has been shown that this may 
be explained in terms of flux capture due to the 
sample boundaries ^B]- Here, Figs. and El sug- 
gest that this effect may persist even for macro- 
scopic homogeneous samples. 
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Figure 4. Vortex configurations in a square (a- 
d) and for a rectangle with L = 8W (e-g), both 
considering W — 0.2A = 20£. The points in the 
magnetization curve corresponding to these con- 
figurations are indicated in Fig. [3] 



6. Conclusions 

In conclusion, we have presented a theoreti- 
cal study of the surface barrier in mesoscopic 
superconducting samples and its effect on the 
formation of metastable vortex structures. Us- 
ing a Langevin dynamics simulation, magneti- 
zation loops of mesoscopic rectangles of several 
aspect ratios were calculated. All these loops 
are hysteretic and present series of magnetization 
jumps due to penetration or expulsion of a vortex. 
An additional modulation of the magnetization 
curve is caused by the formation of vortex chain 
states in the more elongated rectangles. Due to 
the surface barrier, the vortex configurations are 
metastable in both branches of the magnetization 
loop. The capture of vortices in the downward 
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branch leads to paramagnetic response, even for 
the largest samples studied. 
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